Mesoscopic Kondo screening effect in a single-electron transistor embedded in a 

metallic ring 
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We study the Kondo screening effect generated by a single-electron transistor or quantum dot 
embedded in a small metallic ring. When the ring circumference L becomes comparable to the 
fundamental length scale £jj = Hvf/T^ associated with the bulk Kondo temperature, the Kondo 
resonance is strongly affected, depending on the total number of electrons (modulo 4) and magnetic 
flux threading the ring. The resulting Kondo-assisted persistent currents are also calculated in both 
Kondo and mixed valence regimes, and the maximum values are found in the crossover region. 



PACS numbers:72.f5.Qm, 73.23.Ra, 73.23. Hk 

The effect of magnetic impurities on metals — the 
Kondo ..problem has been studied for nearly half a 
centuryEI, and attracts continued interest till now in 
attempts to understand heavy fermion materials and 
high-T c superconductors. Recent progress in the nano- 
fabrication technique of electronic devices has enabled 
investigation of the Kondo effect by means of single- 
electron transistors, £SETs) or quantum dots (QDs) in 
a controllable way0i3. When it contains an odd num- 
ber of localized electrons, a QD may be described by 
the Anderson impurity model, exhibiting the Kondo ef- 
fect characterized by tbe,appearance of a sharp resonance 
near the Fermi energvEra, and it provides a new channel 
for electric current flowing through the QD. 

To investigate the quantum phase interference, a de- 
vice setup of a QD embedded— in a mesoscopic metal- 
lic ring has been proposeo£3~cfl. In such a geometry 
the screening cloud is trapped in the ring, and one can 
measure the persistent current (PC) in the ring induced 
by the external magnetic flux to understand the coher- 
ent transport through the dot without attaching leads 
to it. pBased on perturbation calculations, Affleck and 
SimoncJ recently pointed out that in the Kondo regime, 
the PC is a universal scaling function of where 
= hvp/Tft is the characteristic length scale associ- 
ated with the bulk Kondo temperature T^, and crosses 
over from the value for a perfect ring with no QD in 
the limit ^/I C 1 to a vanishing small magnituderfcr 
^/L> 1, in contrast to the variational calculationsEJ. 

In this Letter, we would like to further clarify the un- 
derlying physics within a model of an Anderson magnetic 
impurity embedded in a mesoscopic metallic ring. Ex- 
plicitly considering the discrete spectrum of conduction 
electrons in the ring using the slave boson mean field 
(SBMF) approach, we show that the Kondo resonance 
is strongly dependent on the ratio of ring circumference 
L to the screening cloud size £^ (L = N, as the lattice 
spacing is assumed to be unity), the total number of elec- 
trons ( modulo 4), and the magnetic flux threading the 



ring. The resulting Kondo-assisted PCs are calculated in 
both Kondo and mixed-valence regimes. Our numerical 
results are in good agreement with those of perturbation 
calculations in the limits of •< 1 and ^/£> lli-i 

The mesoscopic Kondo effect was also studied in ultra- 
small metallic grainsE3. 

SBMF approach. — A QD embedded in a mesoscopic 
Aharonov-Bohm (AB) ring can be described by a one- 
dimensional tight-binding Hamiltonian with N lattice 
sites including an Anderson impurity at site "0" . The 
model Hamiltonian is given by: 
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Here (ti?) denotes the hopping between the QD and 
the left (right) neighboring site of the ring. The phase 
factor (j> is defined by 27r<I>/4>o, where $ and <I>o = h/e 
are external magnetic flux and the flux quantum, re- 
spectively. We only consider a single spin-degenerate 
energy level Ed of the QD. Since the Coulomb repul- 
sion U is generally large in experiments, the limit of 
U — ► oo can be taken. A parameter describing the cou- 
pling strength between the QD and the ring is introduced 
r = ttp(sf) |£(£f)| 2 , where p(sf) and tiep) represent the 
ring density of states (DOS) and hopping at the Fermi 
energy sp , respectively. In this paper, only the half- filled 
case (N e = N, where N e is the total number of localized 
and delocalized electrons) is considered. In the thermo- 
dynamic limit, we would have Sf = 0, p(Q) — N/(2irt), 
* 2 (0) = 2(f£ + tl)/N, and T — (t 2 L + t%)/t. 

The strong coupling SBMF theory has been well- 
known to be a good approximation in describing the 
Kondo physics at low temperatures. Here we will ap- 
ply this method to the Anderson impurity model in the 
presence of a magnetic flux, and we believe that it can 
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still yield reliable results. Because of the infinite U limit, 
the impurity operator can be expressed as S a = /J b in 
the slave-boson representation, where the fermion f a and 
the boson b describe the singly occupied electron and hole 
states, respectively. The constraint b^b + f\f a = 1 
has to be imposed. When the mean field (MF) approx- 
imation is made, the boson operators are replaced by a 
c- number 60 , smd the constraint is implemented by in- 
troducing a chemical potential Ao- Therefore, the MF 
Hamiltonian is written as: 
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H m f — —t 



H Y ( C ia c J+l" + h - C ) + Y Zdflfa 



j=l a 



(tLftc la + tBfPc^U + h.C.) + X (b 2 - 1), (2) 



which is essentially a noninteracting system with a set 
of renormalized parameters id = Ed + Ao and = 
botuny Since the conduction electrons themselves do 
not form a complete set, the Fourier transformation can 
not be used to diagonalize H m f. As a result, the ground 
state energy E gs can only be evaluated by exact diago- 
nalization as 



E, 



occ 



r(Ao,& ) + A (&2_i) ; 



(3) 



where the summation over m includes all occupied levels 
of H m f. The values of Ao and 60 can be determined self- 
consistcntly by minimizing E gs with respect to Ao and 
6 . Once the solutions Ag, 6g and E gs = E* s (\Q,bg) are 
obtained, one can calculate the conduction electron DOS 
at site Ni by 

PNi<j{oj) = -^ Im (0 c Ni<J {uj + ij - H m f)~ l J N . a 0^> , 

and the PC: I = — f^p, where |0) represents the 
vacuum state, and 7 is introduced to denote the half- 
maximum width of the spectrum lines in the DOS. Note 
that the DOS of QD is given by pqd{uj) = bopQ<x(u)), and 
the "Kondo correlation energy" Tk for the finite-size sys- 
tems can be defined by the energy gain due to coupling 
between the QD and the ring: 



Tk - £d - £f - {Eg S - E 



(4) 



where E® s denotes the ground state energy of the nonin- 
teracting Hamiltonian Ho with electron number N e and 
an open boundary condition, and sp is the energy of the 
highest occupied level. Tk is thus the Kondo tempera- 
ture in the finite size system, and it will approach T% in 
the thermodynamical limit [L — > 00). In this sense, we 
will refer to Tk as a mesoscopic Kondo temperature. 

The Anderson impurity model in QD realizations also 
exhibits distinct behaviors in different regimeso. Here 
we will mainly focus on the Kondo regime, namely, the 
highest occupied energy level of the QD has a single 



occupancy. As a typical case, we choose Ed — —0.7t, 
ti, = tft = 0.35i, and T — 0.245t, corresponding to 
T K = 0.0105i or = 191 in the thermodynamic limit. 

DOS of QD. — In the absence of the external mag- 
netic flux, the DOS of QD pqd{u) is obtained by numer- 
ically solving the SBMF self-consistent equations for var- 
ious £,k/L ano - are summarized in Fig. 1, where the half 
maximum width of the spectrum lines has been chosen 
as 7 = 0.5T^-, and different choices of 7 will not change 
the results much. For (J-Ci, the shape of the Kondo 
resonance is indistinguishable from the bulk case. When 
L is decreased down to £ K , however, it splits up into a 
set of subpeaks, in a way which depends strongly on the 
total numbers of lattice sites N (modulo 4). The case 
of N = An + 2 shows two main peaks with nearly the 
same weight around £p, whereas the others only have 
a single main peak at slightly different positions. De- 
spite developing subpeaks, the Kondo resonance retains 
its basic feature, i.e., the profile of the DOS around £p 
is preserved within a few T^- . In other words, the Kondo 
correlations induced by spin-flip processes between the 
singly occupied electron on the QD and the conduction 
electrons persist even in a very small system. 

Features for the ^ K ~ L case are mainly determined by 
the highest occupied and lowest unoccupied single parti- 
cle energy levels of Ti m f- For N — An, when the localized 
electron on QD binds a delocalized electron at the energy 
range |e| < T to form a Kondo singlet near ep — 0, the 
energy difference between the highest occupied and low- 
est unoccupied levels is very small, and a single peak is 
dominant in the DOS. In contrast, for N — An + 2, there 
is a large difference in the resulting highest occupied and 
lowest unoccupied energy levels, so the Kondo resonance 
peaks appear both below and above ef with almost the 
same weight. On the other hand, for the odd parity cases 
of N = An ± 1, the highest energy levels are singly oc- 
cupied below and above Ef , respectively, and a single 
Kondo resonance thus appears below and above Ef sep- 
arately. This feature may have significant consequences. 

The DOS of the QD is also strongly affected by the ex- 
ternal magnetic flux. This is illustrated in Fig. 2, which 
shows the 4> dependence of pqd(uj) for ^ K /L = 1.0. A 
striking feature is that the DOS with N sites and zero 
flux is exactly equal to the one with N + 2 sites and a 
half flux quantum. Such a feature is related to the nature 
of the strong coupling fixed point for the model, where 
one electron is trapped in the symmetric sites around the 
impurity and other delocalized conduction electrons gain 
a resonant phase shift from the scattering off the Kondo 
singlet. It is therefore sufficient to consider only the case 
of N — An and An + 1. For even parity of N = An, the 
single main peak in the DOS gradually splits into two 
symmetric peaks as <f> increases to tt, and they recombine 
again into a single peak as cf> is further increased to 2tt. 
For odd parity of TV = An + 1, the situation is quite 
different. A single main peak in DOS persists asym- 
metrically in the whole range of 4>, and its position is 
shifted with a period of 2tt. However, in this odd par- 
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ity cases there is an important symmetry in the DOS of 
QD: p QD (-uj, 4>) = pq D (u;,it + eft), leading to the period 
halving in the PCs. 

Kondo screening cloud. — In order to determine the 
size of the Kondo screening cloud, Fig. 3a illustrates the 
DOS of the QD and selected lattice sites for a very large 
ring £° K /L = 0.1. At w » 0, the DOS of lattice sites 
Ni near the QD is strongly affected by the Kondo reso- 
nance, and clearly exhibits a parity-dependence, i.e., it 
is enhanced for Ni — even and suppressed for TV, = odd. 
However, such an influence does weaken significantly once 
the lattice site is away from the QD. This is illustrated in 
the change of DOS at sf, as shown in Fig. 3b, which is 
obtained by substracting the background for the nonin- 
teracting part Hq. We observe that ApN ia (oj = 0) drops 
exponentially with increasing Ni as expected. The decay 
length £x deduced from numerical fitting is roughly 190, 
which is close to the Kondo correlation length scale in 
the thermodynamic limit £^ = 191. Note that this is 
not an accurate procedure to determine this characteris- 
tic length scale. 

PC. — The most convenient way to detect the above 
Kondo screening cloud in experiments is to measure the 
PC in a ring induced by the magnetic flux. Such measure- 
ments have been performed recently on micron size rings 
not containing a^DtLj and slightly different experimental 
setup with a QDu. In Fig. 4 we show the PC vs magnetic 
flux for £k/L = 0.3 (dashed line) and 1.5 (solid line). A 
striking feature can be clearly seen I N+2 (<fi) = I n (<P + tt) 
due to the fact that the structures of DOS of QD have an 
exact 7T shift from the system with TV sites to that with 
N + 2 sites. It is obvious that the behavior of PC differs 
greatly in the limit ^/L < 1 from that for £,%/L ~ 1. 
In the former case the system resembles an idea metallic 
ring, implying a high electron transmission through the 
QD. When L becomes comparable to f^-, the PC does 
weaken systematically and their profiles are close to si- 
nusoidal. Moreover, the PC for odd parity appears to be 
much smaller than that for even parity, which also indi- 
cates a weakening of the Kondo screening effect and the 
QD simply acts as an incoherent scatter. Note that these 
features are in good .agreement with the recent pertur- 
bation theory resultsEj, but are partlv in contradiction 
with those of variational calculationst3. 

In order to fully display the dependence on the ring 
circumference, the PCs and Kondo correlation energies 
Tfc are calculated as functions of £,%/L in the cases of 
= 7r/4 and 7r/2, displayed in Fig. 5. We find the PCs in- 
terpolate smoothly between the limits of small and large 
£,%/L. Both I/Iq and Tk/T^ curves follows universal 
scaling functions of the ratio £^ K /L, depending on the 
total number of electrons (modulo 4). It is rather sur- 
prising that the mesoscopic Kondo temperatures Tj< in 
small rings are much higher than the bulk value T^, im- 
plying that the corresponding Kondo correlation lengths 
are much smaller than the bulk value £^ K . 

Finally, we study the properties of the PCs for different 
occupations of the highest occupied energy level of QD by 



tuning the dot energy level Sd- In Fig. 6 the PC is shown 
as a function of Sd for various £%/L at cj> = 7r/4. For dif- 
ferent parities, the PCs exhibit an asymmetric structure 
and their maxima appear near the crossover boundary 
of the Kondo and mixed- valence regimes. In the Kondo 
regime, the PCs are suppressed as L decreases, indicating 
the effects of weakening in the Kondo-assisted tunnelings, 
while in the mixed-valence regime, however, the PCs are 
almost unchanged as L varies. These distinct features 
suggest that the tunneling mechanisms in the Kondo and 
mixed valence regimes are quite different. 

In conclusion, we have calculated the mesoscopic 
Kondo screening effect in a single-electron transistor em- 
bedded in a metallic ring. The Kondo resonance is 
found to be strongly dependent on the dimcnsionless ra- 
tio £j(/L, the total number of electrons (modulo 4), and 
the external flux threading the ring. The screening effects 
on the PCs are expected to be observable in experiments. 
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Figures Captions 

Fig.l. pqd (u>) for various values of Q^/L in the ab- 
sence of a magnetic flux. The number of lattice sites 
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N = L corresponds to (a) 4n, (b) 4n + 1, (c) An + 2 and 
(d) 4n + 3, respectively. Curves are offset by 1.5 units 
each, and the half-maximum width of the spectrum lines 
7 = 0.5T° . 

Fig. 2. pqd (uj) of ^k/L = 1.0 for different external 
flux. In each subplot, from the bottom up, <f> = ~ 2ir 
is increased at steps of 0.2tt. Curves are offset by 2 units 
each. 

Fig. 3. (a) pq D (uj) and selected lattice sites PN t (u) for 
£k/L = 0.1. The dashed lines denote the corresponding 
DOS for Hq, and curves labeled by Ni = 1,2,191, and 
192 are offset by 0.1,0.6,0.8, and 1.1 units, respectively. 
The curve for pqd (oj) is reduced by five times, (b) The 
change of DOS ApN irT (uj — 0) is calculated by substract- 
ing the background for Hq as a function of the lattice 
site Ni . The solid line fits the exponential decay and a 
characteristic length is obtained £k ~ 190. 

Fig. 4. Persistent currents versus flux with the total 
number of lattice sites An+i (i = 0, 1, 2, 3) for (,%/L = 0.3 
(dashed lines) and for £^ K jL = 1.5 (solid lines). Iq de- 
notes the corresponding PC of an ideal metallic ring. 

Fig. 5. Persistent currents and Kondo temperature Tk 
as functions of £%/L for 4> = 7r/4 and tt/2. Both I/I 
and Tk/Tk curves follow universal scaling functions of 
the ratio S^ K /L. 

Fig. 6. Occupation number of the highest occupied en- 
ergy level of QD (solid lines) and PCs I / Iq. The 
dashed, dotted, dash-dotted, and dash-dot-dotted lines 
show the data for = 0.3, 0.6, 1.2, and 2.4, respec- 

tively, as functions of Ed for different parities at (f> = 7r/4. 
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